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Abstract  A  global  high-resolution  ocean  circulation  model  forced  by  atmospheric  fields  and  the  M2  tidal 
constituent  is  used  to  explore  plausible  scenarios  for  the  damping  of  low-mode  internal  tides.  The  plausibil¬ 
ity  of  different  damping  scenarios  is  tested  by  comparing  the  modeled  barotropic  tides  with  TPX08,  a 
highly  accurate  satellite-altimetry-constrained  tide  model,  and  by  comparing  the  modeled  coherent  baro- 
clinic  tide  amplitudes  against  along-track  altimetry.  Five  scenarios  are  tested:  (1)  a  topographic  internal 
wave  drag,  argued  here  to  represent  the  breaking  of  unresolved  high  vertical  modes,  applied  to  the  bottom 
flow  (default  configuration),  (2)  a  wave  drag  applied  to  the  barotropic  flow,  (3)  absence  of  wave  drag,  (4)  a 
substantial  increase  in  quadratic  bottom  friction  along  the  continental  shelves  (with  wave  drag  turned  off), 
and  (5)  application  of  wave  drag  to  the  barotropic  flow  at  the  same  time  that  quadratic  bottom  friction  is 
substantially  increased  along  the  shelves.  Of  the  scenarios  tested  here,  the  default  configuration  (1)  yields 
the  most  accurate  tides.  In  all  other  scenarios  (2-5),  the  lack  of  damping  on  open  ocean  baroclinic  motions 
yields  baroclinic  tides  that  are  too  energetic  and  travel  too  far  from  their  sources,  despite  the  presence  of  a 
vigorous  mesoscale  eddy  field  which  can  scatter  and  decohere  internal  tides  in  the  model.  The  barotropic 
tides  are  also  less  accurate  in  the  absence  of  an  open  ocean  damping  on  barotropic  motions,  that  is,  in  sce¬ 
narios  (3)  and  (4).  The  results  presented  here  suggest  that  low-mode  internal  tides  experience  substantial 
damping  in  the  open  ocean. 
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1.  Introduction 

A  long-standing  problem  in  physical  oceanography  is  the  determination  of  locations  where  tides,  including 
low-mode  internal  tides,  dissipate  their  energy.  An  answer  to  this  problem  will  help  constrain  the  spatial  var¬ 
iability  in  ocean  mixing  patterns  and  consequently  lead  to  more  physically  based  mixing  parameterizations 
for  global  ocean  circulation  models.  It  has  been  argued  that  mixing  driven  by  breaking  internal  tides  and 
other  internal  waves  exerts  a  strong  influence  on  the  stratification  and  general  circulation  of  the  ocean, 
with  consequent  implications  for  climate  variability  [ Munk ,  1966;  Munk  and  Wunsch,  1998;  Wunsch  and 
Ferrari,  2004]. 

Using  barotropic  tidal  models  constrained  by  TOPEX/Poseidon  satellite  altimeter  data,  Egbert  and  Ray  [2000, 
2001]  clearly  demonstrated  that  a  significant  percentage  (25-30%)  of  barotropic  tidal  energy  is  lost  in  open 
ocean  regions  of  rough  topography.  The  energy  loss  in  the  Egbert  and  Ray  studies  represents  the  net  effect 
of  mechanisms  that  siphon  energy  from  the  barotropic  tide  including  conversion  to  baroclinic  tides  as  well 
as  loss  to  quadratic  bottom  boundary  layer  drag  (hereafter  referred  to  as  "quadratic  bottom  friction"). 
Because  internal  tides  were  not  directly  simulated  in  the  Egbert  and  Ray  studies,  they  could  not  discriminate 
between  barotropic  tidal  energy  loss  due  to  local  dissipation  and  due  to  generation  and  radiation  of  internal 
tides.  Using  along-track  data  from  the  TOPEX/Poseidon  satellite-altimetry  mission,  Ray  and  Mitchum  [1996, 
1997]  revealed  the  surface  manifestations  of  the  coherent  low-mode  internal  tides  that  propagate  up  to 
1 000  km  from  sources  such  as  the  Hawai'ian  Ridge  [also  see  Alford  and  Zhao,  2007;  Zhao  et  al.,  2010].  The 
long-range  propagation  of  internal  tides  from  sources  such  as  Hawai'i  has  also  been  revealed  in  an  analysis 
of  acoustic  tomography  [Dushaw  et  al.,  1995]  and  an  analysis  of  tide  gauge  data  supplemented  with 
inverted  echo  sounders  [Mitchum  and  Chiswell,  2000].  Where  and  how  these  low-mode  internal  tides  are 
damped  is  a  topic  of  great  interest  in  the  community. 
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Direct  simulation  of  baroclinic  tides  has  been  done  on  regional  scales  [e.g.,  Carter  et  al.,  2012,  and  references 
therein]  and  basin  scales  [e.g.,  Niwa  and  Hibiya,  2001],  as  well  as,  more  recently,  global  scales.  Using  three- 
dimensional  ocean  models  forced  by  only  the  astronomical  tidal  potential,  Arbic  et  al.  [2004]  and  Simmons 
et  al.  [2004]  directly  simulated  the  global  baroclinic  tides,  further  revealing  the  rich  structure  of  the  internal 
tide  climate  in  the  ocean.  Arbic  et  al.  [2010]  later  showed  that  this  rich  structure  can  still  be  captured  in  a 
model  that  is  simultaneously  forced  by  both  the  astronomical  tidal  potential  and  atmospheric  fields  [see 
also  Arbic  et  al.,  2012;  Muller  et  al.,  2012;  Waterhouse  et  al.,  2014]  (C.  B.  Rocha,  Mesoscale  to  submesoscale 
wavenumber  spectra  in  Drake  Passage,  submitted  to  Journal  of  Physical  Oceanography,  2015,  based  on  sim¬ 
ulations  by  D.  Menemenlis  (personal  communication,  2014,  2015)).  Because  high-resolution  global  models 
with  simultaneous  atmospheric  and  tidal  forcing  are  expensive  and  relatively  new,  no  comprehensive  inves¬ 
tigation  has  been  done  thus  far  of  internal  tide  damping  in  such  models. 

Several  mechanisms  have  been  identified  in  the  literature  as  potential  dissipators  of  the  internal  tides.  After 
generation,  the  high-vertical-mode  internal  tides  may  break  locally  [Klymak  et  al.,  2006,  2008;  Legg  and  Kly- 
mak,  2008;  Buijsman  et  al.,  2012]  while  the  remaining  low-mode  tides  propagate  away  from  their  topo¬ 
graphic  sources.  Propagating  low-mode  tides  may  dissipate  via  interactions  with  distant  topographies 
[Johnston  and  Merrifield,  2003;  Mathur  et  al.,  2014]  or  through  interactions  with  mean  flows  and  eddies  [Sf. 
Laurent  and  Garrett,  2002;  Rainville  and  Pinkel,  2006;  Dunphy  and  Lamb,  2014;  Kerry  et  al.,  2014].  They  may 
also  lose  energy  through  upper-ocean  wave-wave  interactions,  for  instance  through  the  parametric  subhar¬ 
monic  instability  (PSI)  mechanism  [McComas  and  Bretherton,  1 977;  Muller  et  al.,  1 986;  Staquet  and  Sommeria, 
2002;  MacKinnon  and  Winters,  2005;  Simmons,  2008;  Sun  and  Pinkel,  2012,  2013].  Finally,  low-mode  internal 
tides  may  propagate  onto  continental  shelves  and  subsequently  dissipate,  thus  mixing  the  coastal  ocean 
[Nash  et  al.,  2004;  Martini  et  al.,  2011;  Kelly  et  al.,  2013;  Waterhouse  et  al.,  2014].  Detailed  discussions  of  the 
above  mechanisms  can  be  found  in  previous  publications  [e.g.,  St.  Laurent  and  Garrett,  2002;  Kunze  and  Lle¬ 
wellyn  Smith,  2004],  The  relative  importance  of  each  of  these  mechanisms  in  the  global  ocean  is  unclear  at 
present. 

Global  ocean  models  have  become  indispensable  in  understanding  large-scale  oceanic  motions,  due  to 
their  widespread  spatial  coverage,  relatively  frequent  output  and  relatively  low  cost  compared  with  obser¬ 
vational  systems.  Nevertheless,  global  models  have  many  shortcomings.  For  instance,  global  models  are  not 
able  to  resolve  the  breaking  of  internal  gravity  waves.  Justifications  for  the  usage  of  a  parameterized  wave 
drag  in  global  barotropic  tide  models  are  expounded  upon  in  several  studies  [e.g.,  Jayne  and  St.  Laurent, 
2001;  Arbic  et  al.,  2004;  Egbert  et  al.,  2004;  Green  and  Nycander,  2013;  Buijsman  et  al.,  2015].  Among  others, 
Arbic  et  al.  [2004,  2010]  argued  that  the  generation  of  low-mode  internal  tides  in  current  state-of-the-art 
global  internal  tide  models  does  not  preclude  the  need  for  an  explicit  damping  of  low-mode  internal  tides 
in  such  models,  which  are  unable  to  resolve  generation  and  breaking  of  high-mode  internal  waves.  In  Arbic 
et  al.  [2004,  2010]  and  in  subsequent  related  papers  [e.g.,  Shriver  et  al.,  2012;  Buijsman  et  al.,  2015],  it  was 
argued  that  a  parameterized  wave  drag  acting  on  the  bottom  flow  can  be  taken  to  represent  the  genera¬ 
tion  and  breaking  of  unresolved  high  modes  by  flow  over  rough  topography.  Simultaneously  Shriver  et  al. 
[2012]  showed  that  the  barotropic  and  low-mode  baroclinic  tides  from  simulations  forced  by  tides  and 
atmospheric  fields  compare  well  to  satellite  altimeter  observations  when  a  parameterized  internal  wave 
drag  is  applied  to  the  bottom  flow.  Note  that  because  the  bottom  flow  has  both  barotropic  and  low-mode 
baroclinic  components,  such  a  drag  damps  low-mode  internal  tides  as  well  as  barotropic  tides. 

We  can  test  some  of  the  aforementioned  low-mode  internal  tide  damping  mechanisms,  by  applying  them 
to  global  internal  tide  simulations  and  comparing  the  resulting  barotropic  and  baroclinic  tides  with 
altimetry-constrained  models  and  observations,  as  was  done  in  Shriver  et  al.  [2012]  for  a  wave  drag  acting 
on  the  bottom  flow.  By  comparing  the  resulting  tides  to  observations,  we  can  determine  which  drag  mech¬ 
anisms  are  plausible  and  which  are  not.  In  the  scenarios  we  are  testing,  we  will 

a.  apply  parameterized  topographic  internal  wave  drag  on  the  bottom  flow  (referred  to  as  the  "default 
configuration"  or  "bottom  wave  drag"), 

b.  apply  wave  drag  on  only  the  barotropic  tides  (referred  to  as  "barotropic  wave  drag"), 

c.  remove  the  wave  drag  entirely, 

d.  remove  wave  drag  and  substantially  increase  quadratic  bottom  friction  along  the  shelves,  and 

e.  apply  wave  drag  on  only  the  barotropic  tides  and  at  the  same  time  substantially  increase  quadratic 
bottom  friction  along  the  shelves. 
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Scenario  (a)  represents  scattering  into  high  modes  above  rough  topography.  We  note  that  we  will  vary  the 
strength  of  the  wave  drag  in  scenario  (a)  but  will  refer  to  all  such  simulations  with  wave  drag  applied  on  the 
bottom  flow  as  "default  configuration"  or  "bottom  wave  drag."  In  scenario  (b),  only  the  barotropic  mode 
experiences  topographic  scattering,  and  again  we  will  vary  the  drag  strength  in  this  scenario.  In  scenario  (c), 
neither  barotropic  nor  low-mode  baroclinic  tides  experience  substantial  damping  in  the  open  ocean.  Sce¬ 
nario  (d)  represents  the  hypothesis  that  most  damping  of  low-mode  tides  occurs  in  shallow  waters.  Scenario 
(e)  is  a  combination  of  (b)  and  (d)  that  prevents  overly  energetic  generation  of  baroclinic  tides  in  the  open 
ocean  while  still  testing  whether  most  damping  of  low-mode  internal  tides  occurs  along  the  shelves.  The 
scenarios  (a)-(e)  are  not  an  exhaustive  list  of  potential  low-mode  internal  tide  damping  mechanisms.  Damp¬ 
ing  by  upper-ocean  wave-wave  interactions  is  not  tested  here  and  will  instead  be  tested  in  a  paper  under 
preparation,  in  which  simplified  global  tidal  models  (with  only  tidal  forcing)  will  be  used.  Furthermore, 
damping  of  low-mode  tides  in  the  actual  ocean  is  likely  to  arise  from  more  than  one  mechanism.  However, 
it  is  too  complicated  and  expensive  in  this  first  study  to  utilize  complex  combinations  of  different  mecha¬ 
nisms  in  our  simulations.  Our  simulations  represent  various  limiting  cases  and  will  likely  point  further  stud¬ 
ies  in  fruitful  directions.  We  further  note  that  because  energy  budget  analyses  of  our  simulations  are 
complex  and  computationally  expensive,  we  do  not  undertake  them  here  and  instead  will  report  on  them 
elsewhere.  A  detailed  energy  budget  analysis  for  one  of  our  previous  HYCOM  simulations  is  in  review  (M.  C. 
Buijsman  et  al.,  Impact  of  internal  wave  drag  on  the  semidiurnal  energy  balance  in  a  global  ocean  circula¬ 
tion  model,  submitted  to  Journal  of  Physical  Oceanography,  2015).  Our  focus  in  the  present  paper  is  on  test¬ 
ing  the  plausibility  of  extreme  damping  scenarios. 

It  should  be  noted  that  in  all  scenarios  examined  in  this  paper,  the  presence  of  atmospheric  forcing  in  a 
high-resolution  model  ensures  a  vigorous  mesoscale  eddy  field  [Hecht  and  Hasumi,  2008;  Thoppil  et  al., 
201 1].  The  mesoscale  eddies  will  decrease  the  coherence  of  the  internal  tides,  or  in  other  words,  will  render 
a  nonstationary  component  of  the  internal  tide  field  [e.g.,  Shriver  et  al.,  2014].  Because  along-track  altimetry 
only  detects  the  coherent  component  of  baroclinic  tides,  it  has  been  argued  that  the  general  weakness  of 
low-mode  tides  far  from  their  sources  seen  in  global  maps  made  from  altimetry  could  reflect  tidal  incoher¬ 
ence  arising  from  scattering  by  eddies,  rather  than  damping.  Because  eddy  scattering  is  present  in  our  sim¬ 
ulations,  the  presence  of  baroclinic  tides  in  our  model  that  are  stronger  than  or  propagate  farther  than  tides 
in  observations  will  be  taken  as  evidence  of  under-damping  in  the  model,  rather  than  a  signature  of  an 
under-represented  scattering,  as  could  be  the  case  in  a  tide-only  model. 

To  our  knowledge,  global  baroclinic  tide  simulations  representing  cases  (b),  (d),  and  (e)  have  not  been 
reported  before.  Several  results  [ Arbic  et  al.,  201 0,  201 2;  Richman  et  al.,  2012;  Shriver  et  al.,  2012,  2014;  Timko 
et  al.,  2012,  2013;  Muller  et  al.,  2015]  have  been  published  on  case  (a)  using  simulations  of  the  Hybrid  Coor¬ 
dinate  Ocean  Model  (HYCOM)  [Bleck,  2002;  Chassignet  et  al.,  2003;  Halliwell,  2004].  The  Arbic  et  al.  [2010]  sim¬ 
ulation  was  preliminary  and  has  not  been  used  in  any  of  our  subsequent  papers.  The  case  (a)  scenario  used 
in  the  other  HYCOM  papers  mentioned  above  was  run  for  a  long  time  and  will  prove  useful  in  this  paper, 
which  focuses  most  of  its  attention  on  runs  performed  for  the  first  time  here.  Global  baroclinic  tide  simula¬ 
tions  in  some  previous  studies  [ Simmons  et  al.,  2004;  Simmons,  2008;  Muller  et  al.,  2012;  Waterhouse  et  al., 
2014]  were  performed  without  utilizing  a  parameterized  topographic  wave  drag  (referred  to  as  case  c 
above).  To  obtain  accurate  barotropic  tides  in  their  basin  and  global  models,  Niwa  and  Hibiya  [2001,  201 1] 
employed  a  linear  damping  scheme,  representing  energy  loss  via  nonlinear  interactions  with  the  back¬ 
ground  internal  waves  and  mesoscale  eddies,  to  account  for  the  decay  of  internal  tides  as  they  propagate 
from  their  sources. 

In  the  next  section,  we  present  the  HYCOM  simulations  utilized  in  this  paper  and  the  diagnostics  we  will 
employ  to  compare  the  model  output  with  observations.  Section  3  presents  the  results  of  the  simulated 
barotropic  and  baroclinic  tides  and  their  comparison  with  observations  under  the  different  drag  scenarios. 
Section  4  presents  the  summary  and  conclusions. 

2.  Model  and  Diagnostics 

2.1.  The  Model 

This  study  utilizes  simulations  of  HYCOM,  which  is  used  by  the  United  States  Navy  as  an  operational  model. 
The  HYCOM  simulations  presented  here  are  run  in  forward  (non-data-assimilative)  mode  and  utilize  41 
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Table  1.  Description  of  the  Different  Simulations  and  the  Strength  of  the  Applied 
Wave  Drag  and  Quadratic  Bottom  Friction3 

Drag  Scale 


Simulation 

Description 

Factor 

E053 

No-wave  drag 

0.0 

E054 

Increased  bottom  friction  (X25)  along  shelves 

0.0 

E055 

Increased  bottom  friction  (XI 00)  along  shelves 

0.0 

E072 

Wave  drag  on  bottom  flow 

0.2 

E073 

Wave  drag  on  barotropic  flow 

0.2 

E051 

Wave  drag  on  bottom  flow 

0.5 

E052 

Wave  drag  on  barotropic  flow 

0.5 

E074 

Wave  drag  on  barotropic  flow  and  increased 
bottom  friction  (XI 00)  along  shelves 

0.5 

E056 

Wave  drag  on  bottom  flow 

0.7 

E057 

Wave  drag  on  barotropic  flow 

0.7 

E058 

Wave  drag  on  bottom  flow 

1.0 

E059 

Wave  drag  on  barotropic  flow 

1.0 

hybrid  vertical  layers  on  a  global  tri¬ 
pole  grid  with  a  nominal  horizontal 
resolution,  at  the  equator,  of  1/ 
12.5°.  The  simulations  are  forced  by 
the  M2  tide,  the  largest  tidal  constit¬ 
uent  in  the  ocean,  and  by  hourly 
atmospheric  fields  from  the  Navy 
Global  Environmental  Model  (NAV- 
GEM)  [Hogan  et  al.,  2014].  The  initial 
state  is  from  a  case  with  HYCOM 
coupled  to  the  Los  Alamos  sea  ice 
model  (CICE)  [Hunke  and  Lipscomb, 
2008],  initialized  from  a  January  cli¬ 
matology,  spun  up  with  climatologi¬ 
cal  atmospheric  forcing  for  7  model 
years  and  then  run  from  2003  to 
June  201 1  using  Navy  Operational 
Global  Atmospheric  Prediction  Sys¬ 
tem  (NOGAPS)  [Rosmond  et  al.,  2002]  atmospheric  forcing.  NAVGEM  replaced  NOGAPS  as  the  US  Navy's 
operational  global  numerical  weather  prediction  system  in  2013.  Our  simulations  began  to  use  NAVGEM  in 
June  201 1  because  we  have  hourly  NAVGEM  fields  for  that  time  period.  The  June  201 1  ocean  from  HYCOM 
was  first  extended  under  the  Antarctic  ice  shelf,  to  improve  tidal  accuracy,  and  CICE  was  replaced  by  HYCOM's 
built-in  thermodynamic  sea  ice  model  to  save  computer  time.  HYCOM  has  several  options  for  tidal  Self  Attrac¬ 
tion  and  Loading  (SAL)  [ Hendershott ,  1972].  Here  we  use  the  observed  SAL,  calculated  from  the  data- 
assimilative  TPX08-atlas  [Egbert  and  Erofeeva,  2002,  updated],  for  all  cases.  Most  of  the  simulations  here 
are  run  from  July  201 1  to  March  2012.  Hourly  sea  surface  height  (SSH)  data  were  saved  for  the  last  3  months 
(January-March  2012)  of  each  simulation.  Four  simulations  were  run  for  a  longer  period  of  10  months  to 
investigate  the  dependence  of  the  results  upon  the  length  of  time  series  as  will  be  discussed  in  section  3.2.2. 


aln  E054,  E055,  and  E074,  the  bottom  friction  (Q)  along  the  shelves  (depths 
<2500  m)  is  increased  by  about  25  and  100  times,  respectively,  above  the  value  in 
the  deep  ocean  (Q  =  0.0025)  (see  text  for  additional  explanation). 


In  this  paper,  we  employ  the  topographic  wave  drag  scheme  by  Jayne  and  St.  Laurent  [2001].  The  drag  was 
tuned  in  baroclinic  HYCOM  simulations  with  different  scale  factors  showing  that  a  scale  factor  of  0.5  is  opti¬ 
mum.  Other  scale  factors  will  also  be  utilized  in  our  simulations.  Our  previous  HYCOM  simulations  [e.g., 
Shriver  et  al.,  2012]  used  the  scheme  by  Garner  [2005].  Here  we  turn  to  the  Jayne  and  St.  Laurent  [2001] 
scheme  because  of  its  simplicity  and  smoother  nature  [Green  and  Nycander,  2013;  Buijsman  et  al.,  2015] 
which  makes  it  easier  to  use  in  baroclinic  models  and  less  prone  to  numerical  instabilities.  We  performed 
two  additional  simulations  with  the  Garner  scheme  showing  that  it  also  gives  good  barotropic  and  baro¬ 
clinic  tides.  Later,  we  will  also  turn  to  an  older  simulation  (designated  internally  as  HYCOM  18.5)  that 
employed  the  Garner  [2005]  drag  scheme,  and  that  has  been  used  in  many  of  our  HYCOM  tides  papers  to 
date  [e.g.,  Shriver  et  al.,  2012,  among  others].  The  older  simulation  was  run  out  for  5  years  and  will  be  useful 
for  our  analysis  of  the  impact  of  record  durations  at  the  end  of  the  results  section.  A  brief  discussion  of  the 
wave  drag  schemes  by  Jayne  and  St.  Laurent  [2001]  and  Garner  [2005]  is  given  below.  Thorough  discussions 
on  topographic  wave  drag  and  quadratic  bottom  friction  and  their  appearance  in  the  momentum  equations 
can  be  found  in  previous  publications  [e.g.,  Arbic  et  al.,  2004,  2010;  Buijsman  et  al.,  2015].  In  particular,  Arbic 
et  al.  [2010]  discuss  the  subtleties  involved  in  applying  wave  drag  to  only  the  tidal  flow  in  HYCOM  simula¬ 
tions  that  simultaneously  resolve  both  tidal  and  nontidal  motions.  Table  1  displays  the  12  different  simula¬ 
tions  conducted  for  this  paper  with  additional  details  to  be  explained  below. 


2.2.  Topographic  Wave  Drag  Schemes 

The  topographic  wave  drag  scheme  by  Jayne  and  St.  Laurent  [2001]  is  represented  by  a  linear  drag  coefficient 
that  varies  spatially  based  upon  the  topography  [also  see  Green  and  Nycander,  2013;  Buijsman  et  al.,  2015],  viz: 

CjSL=oc^H2Nb,  (1) 

where  H  is  the  bottom  roughness,  Nb  is  the  buoyancy  frequency  at  the  bottom,  and  L  is  the  wave  length  of 
the  topography.  L  is  set  to  10  km  and  a  =  1  in  Jayne  and  St.  Laurent  [2001].  This  bottom  wave  drag 
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Figure  1.  Map  of  variable  bottom  friction  around  the  Gulf  of  Mexico,  using  (a)  the  nominal  choice  of  Zo  =  10  mm  and  (b)  Zo  =  1 .13  m  (resulting  in  quadratic  drag  coefficient  that  can  be 
about  100  times  larger  than  the  deep  ocean  value  of  0.0025).  Note  the  different  scales  used  in  the  two  subplots. 


coefficient  acts  upon  flow  in  the  bottom  500  m,  which  is  taken  as  a  typical  vertical  scale  length  for  internal 
tides  following  discussions  in  St.  Laurent  and  Garrett  [2002].  In  HYCOM,  a  is  used  as  a  wave  drag  scale  factor 
(or  a  tuning  parameter)  to  tune  the  model  for  accurate  barotropic  and  baroclinic  tides  with  respect  to  obser¬ 
vations.  Different  values  of  a  are  utilized  in  the  simulations  presented  in  this  paper. 

The  scheme  by  Garner  [2005]  is  based  on  an  exact  analytical  solution  for  drag  in  the  limit  of  arbitrary 
small-amplitude  topography.  The  scheme  is  supplemented  by  scaling  arguments  to  account  for  nonlinear 
blocking  effects.  As  with  the  Jayne  and  St.  Laurent  [2001]  scheme,  the  Garner  [2005]  scheme  takes  the 
buoyancy  frequency  at  the  bottom,  and  the  seafloor  topographic  field,  as  inputs.  For  simplicity,  the  origi¬ 
nal  tensor  scheme  in  Garner  [2005]  is  reduced  to  a  scalar  form  in  our  FIYCOM  simulations.  As  with  the 
Jayne  and  St.  Laurent  [2001]  scheme,  the  Garner  [2005]  scheme  is  applied  over  the  bottom  500  m  of  the 
flow.  Arbic  et  al.  [2010]  provide  more  details  of  the  Garner  [2005]  scheme  and  its  implementation  in 
HYCOM. 

The  strength  of  the  wave  drag  can  be  measured  by  computing  the  spatially  averaged  decay  rate.  In  the 
case  of  the  scheme  by  Jayne  and  St.  Laurent  [2001  ],  this  is  given  by  [Arbic  et  al.,  2004] 


_JJ(<WH)dA 

^  JJcM  ' 


(2) 


where  H  is  the  resting  water  depth  and  dA  is  an  element  of  area.  We  computed  the  spatially  averaged  decay 
rate  in  waters  deeper  than  1000  m  for  wave  drag  scale  factors  of  0.5  and  1.0  and  obtained  (3.8days)_1  and 
(1 .9days)~\  respectively.  A  similar  calculation  for  our  older  simulation  using  the  Garner  [2005]  scheme  with 
a  drag  scale  factor  of  12  gives  (1.3days)~\  We  note  that  these  numbers  are  comparable  to  the 
(~  2.0days)~1  decay  rate  computed  by  Arbic  et  al.  [2004]  for  their  barotropic  tide  simulations. 


2.3.  Bottom  Friction  Along  the  Shelves 

In  the  typical  HYCOM  configuration  [see  Buijsman  et  al.,  2015],  the  quadratic  drag  coefficient  is  0.0025  in 
deep  water  and  is  a  function  of  depth  [Schlichting,  1968;  Oey,  2006] 


Q  = 


K 

log  (0.5 H/zo) 


(3) 


in  shallow  water,  where  k=0.4  is  the  von  Karman  constant,  z0  is  the  roughness  parameter,  and  H  is  the  rest¬ 
ing  water  depth.  The  nominal  choice  of  z0  =  1 0  mm  gives  values  of  Q  that  are  about  twice  as  high  as  0.0025 
for  depths  shallower  than  about  60  m. 
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To  substantially  increase  the  bottom  friction  along  the  shelves  for  cases  (d)  and  (e)  described  earlier,  we 
restrict  Q=0.0025  for  all  locations  having  seafloor  depths  greater  than  2500  m.  We  then  choose  two  differ¬ 
ent  values  of  z0  (0.5  and  1.13  m)  to  determine  the  values  of  Q  along  the  shelves  using  equation  (3).  The  first 
choice  z0=0.5  m  results  in  Cd  that  can  be  about  25  times  higher  than  0.0025  around  the  shelves  (see  Figure 
1  for  an  example  plot  around  the  Gulf  of  Mexico)  while  the  second  choice  results  in  Cd  that  can  be  about 
100  times  above  0.0025.  The  first  choice  is  utilized  in  simulation  E054  of  Table  1,  while  the  second  choice  is 
utilized  in  simulations  E055  and  E074. 

2.4.  Diagnostics 

We  evaluate  the  accuracy  of  tides  in  the  different  simulations  by  computing  the  signal  size  (S)  and  the  root- 
mean-square  error  (RMSE)  from  time  series  of  SSH.  The  barotropic  tides  are  validated  against  the  TPX08- 
atlas  sea  surface  elevations  while  the  baroclinic  tides  are  compared  to  almost  17  years  of  along-track 
TOPEX/Poseidon  and  Jason  satellite  altimeter  data  [Ray  and  Mitchum,  1996;  Ray  and  Byrne,  2010;  Shriver 
et  ai,  2012;  R.  Ray,  personal  communication,  2011].  Details  of  the  analysis  approach  used  to  extract  baro¬ 
clinic  tides  from  the  HYCOM  output  and  along-track  altimeter  data  are  outlined  in  Shriver  et  ai.  [2012].  In 
brief,  FIYCOM  output  is  interpolated  to  the  altimeter  tracks.  For  both  FIYCOM  and  altimeter  output,  spatial 
band-pass  filtering  along  the  altimeter  tracks  is  used  to  extract  the  M2  internal  tide  signals  with  wavelengths 
in  the  50-400  km  range.  Only  locations  with  sea  floor  depth  greater  than  1500  m  are  examined.  The  RMSE 
is  defined  as  [Arbic  et  ai,  2004] 


and 

RMSEt  =  ^/  <Dlm(f)  — »/0(f)]2},  (5) 

where  r\  is  sea  surface  elevation,  subscripts  “m"  and  “o"  refer  to  model  and  observations,  respectively,  f  is 
time,  <  >  represents  time  averaging,  and  dA  is  an  element  of  area.  By  writing  the  modeled  and  observed 
elevations  as  r/m(f)=Amcos  (wt-cf>m)  and  r/0(f)  =A0cos  (cot- 4>a),  where  A  and  </>  are  the  amplitude  and 
phase  of  the  signals,  respectively,  RMSEt  greatly  simplifies  to  [e.g.,  Cummins  and  Oey,  1997] 


RMSE,^^(Am-A0)2+AmA0[: -cos  (<£m-</>0)].  (6) 

The  signal  of  the  model  elevations  is  defined  as  [Arbic  et  ai.,  2004] 


c_  /.[.[<  K(t)}2>  dA_  /nxcw 

V  JJdA  v  J7^  ’ 

while  the  signal  of  the  observations  is  similar,  with  r/m  and  Am  replaced  by  i]a  and  Aol  respectively. 

3.  Results 

3.1.  Barotropic  Tides 

Amplitude  and  phase  maps  of  the  M2  surface  elevations  in  TPX08  and  FIYCOM  are  shown  in  Figure  2.  Surface 
elevations  in  FIYCOM  include  the  internal  tide  but  are  dominated  by  the  barotropic  tide.  The  FIYCOM  simula¬ 
tions  displayed  in  Figures  2b-2d  show  qualitatively  similar  spatial  distributions  to  those  in  TPX08.  They  cap¬ 
ture  almost  all  amphidromic  locations  correctly  except  at  the  Antarctic  coast  around  160°W  where  Figures  2c 
and  2d  deviate  from  TPX08.  A  notable  difference  occurs  around  the  Amazon  and  the  east  coast  of  America 
where  FIYCOM  underpredicts  tidal  amplitudes.  FIYCOM  simulations  in  which  the  wave  drag  is  completely 
turned  off  (Figure  2e),  and  in  which  there  is  a  strong  increase  of  quadratic  bottom  friction  along  the  shelves 
but  with  a  turned-off  wave  drag  (e.g.,  Figure  2f),  produce  barotropic  tides  that  are  overly  energetic  relative  to 


ANSONG  ETAL. 


OPEN  OCEAN  INTERNAL  TIDE  DAMPING 


6062 


©AGU  Journal  of  Geophysical  Research:  Oceans 


10.1 002/201 5JC01 0998 


80°N 

60°K 

40°N 

20°N 

0° 

20°S 

40°S 

60°S 

80£s 
80  N 

60°N 

40°N 

20°N 

0° 

20°S 

40°S 

60°S 

80°S 

80°N 

60°N 

40°N 

20°N 

0° 

20°S 

40°S 

60°S 

80°S 


Longitude  (degrees) 


Longitude  (degrees) 


Figure  2.  Amplitude  (colors)  and  phase  (white  contours  separated  by  30°)  of  the  M2  surface  elevation  in  (a)  TPX08-atlas,  a  barotropic  tide  model  constrained  by  satellite  altimetry,  and 
HYCOM  simulations  (b)  E051;  with  wave  drag  (scale  factor  =  0.5)  applied  to  the  bottom  flow,  (c)  E058;  with  wave  drag  (scale  factor  =  1.0)  applied  to  the  bottom  flow,  (d)  E059;  with  wave 
drag  (scale  factor  =  1.0)  applied  to  only  the  barotropic  flow,  (e)  E053;  without  wave  drag,  (f)  E055;  without  wave  drag  but  with  quadratic  bottom  friction  increased  by  about  100  times 
along  the  continental  shelves. 


TPX08.  This  is  especially  visible  in  the  Pacific  and  Indian  oceans,  although  the  amplitudes  in  the  Amazon 
region  appear  more  similar  to  observations  than  they  do  in  the  simulations  with  wave  drag. 

Figure  3  (black  curves)  displays  the  RMSE  in  the  HYCOM  simulations  measured  against  TPX08  showing  that 
the  bottom  wave  drag  simulation  E051  (with  wave  drag  scale  factor  of  0.5)  yields  the  most  accurate  baro¬ 
tropic  tides.  For  any  given  scale  factor  of  the  wave  drag,  we  see  that  more  accurate  barotropic  tides  are 
obtained  if  the  wave  drag  is  applied  to  the  bottom  flow  (bottom  wave  drag;  simulations  E051,  E056,  E058, 
and  E072).  Application  of  the  wave  drag  to  only  the  barotropic  component  of  the  flow  (barotropic  wave 
drag;  simulations  E052,  E057,  E059,  and  E073)  results  in  somewhat  less  accurate  barotropic  tides  compared 
to  their  counterparts  with  the  bottom  wave  drag.  Completely  removing  the  wave  drag  (simulation  E053) 
increases  the  barotropic  SSH  error  by  about  2.0  cm  above  that  of  the  most  accurate  bottom  wave  drag  case 
(simulation  E051).  This  is  because  the  absence  of  a  wave  drag  results  in  overly  energetic  barotropic  tides  as 
seen  in  Figure  2e.  Figure  3  also  shows  that  increasing  the  quadratic  bottom  friction  along  the  shelves,  with¬ 
out  including  a  wave  drag,  yields  the  largest  RMSE  values  (E054  and  E055).  Simulation  E074  is  similar  to 
E052  (both  have  a  barotropic  wave  drag)  except  that  E074  has  large  quadratic  bottom  friction  along  the 
shelves.  The  additional  increase  in  bottom  friction  in  E074  yields  slightly  less  accurate  tides.  Note  that  the 
plots  of  RMSE  in  Figure  3  display  results  for  three  different  months  of  model  output  from  each  simulation. 
We  find  that  there  is  very  little  monthly  variance  in  the  barotropic  tidal  RMSE.  We  performed  a  similar  RMSE 
analysis  on  the  barotropic  tides  using  3  month-long  SSH  data  and  found  that  the  results  were  little  different 
from  the  1  month  analysis  results.  This  contrasts  with  our  findings  in  the  case  of  the  internal  tides,  to  be  dis¬ 
cussed  later. 

The  magenta  curves  in  Figure  3  display  the  signal  S  of  each  simulation  and  that  of  TPX08.  We  find  that 
the  signal  generally  decreases  with  increasing  scale  factor,  because  more  energy  is  removed  from 
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Figure  3.  The  root-mean-square  error  (RMSE;  black  curves)  of  HYCOM  surface 
tidal  elevations  measured  against  TPX08.  The  mean  signal,  5,  in  HYCOM  and 
TPX08  is  depicted  in  the  magenta  curves  (right  axis).  The  scale  factor  of  the  wave 
drag  for  each  simulation  is  shown  on  the  top  axis  and  the  bottom  and  barotropic 
drag  simulations  are  denoted  by  red  and  blue  vertical  dotted  lines,  respectively. 


the  barotropic-baroclinic  system.  As 
expected,  the  simulation  with  the  small¬ 
est  RMSE  (E051)  also  has  an  S  value  that 
lies  close  to  the  TPX08  value.  The  simu¬ 
lations  without  a  wave  drag  (simula¬ 
tions  E053,  E054,  and  E055),  even  those 
with  enhanced  quadratic  bottom  fric¬ 
tion  along  the  shelves,  have  excessively 
large  barotropic  amplitudes. 

A  surprising  outcome,  at  least  to  us,  of 
the  barotropic  tidal  analysis  is  how  small 
the  differences  in  the  errors  in  Figure  3 
are  (<2.2  cm).  As  opposed  to  barotropic 
models  [Arbic  et  at.,  2004;  Jayne  and  St. 
Laurent,  2001;  Buijsman  et  at.,  201 5],  and 
baroclinic  models  with  only  two  layers 
[Arbic  et  at.,  2004],  that  show  a  large 
range  in  errors  depending  on  whether 
wave  drag  is  utilized  or  not,  higher  reso¬ 
lution  baroclinic  models  with  more  verti¬ 
cal  layers  appear  to  show  smaller  ranges. 


3.2.  Baroclinic  Tides 

3.2.1.  Steric  Sea  Surface  Height  Analyses 

Before  we  compare  the  simulated  global  internal  tides  to  along-track  altimeter  data,  we  first  present  the 
internal  tides  from  the  different  scenarios.  For  this  first  look  at  internal  tides,  we  use  the  steric  SSFH  fields. 
The  steric  SSFI  refers  to  changes  in  sea  level  due  to  thermal  expansions  and  salinity  variations  (i.e.,  the  sea 
surface  height  component  having  to  do  with  stratification  effects).  The  steric  SSFI  fields  cover  the  entire 
globe  and  are  direct  model  outputs  without  any  contamination  due  to  spatial  band  passing  and/or  interpo¬ 
lation  of  results  to  satellite  tracks.  The  drawback  of  using  steric  SSFI  fields  is  the  unavailability  of  observatio¬ 
nal  data  for  comparison. 

The  amplitude  of  the  M2  internal  tide  signature  in  the  steric  sea  surface  height,  recovered  from  a  tidal  har¬ 
monic  analysis  of  3  months  of  the  model  steric  fields,  is  displayed  in  Figure  4.  The  beam-like  structures  of 
the  tides  radiating  away  from  their  generation  sources  [e.g.,  Simmons  et  at.,  2004]  are  readily  apparent.  The 
beams  are  generated  in  hot  spot  locations,  which  are  depicted  in  Figure  4  by  the  rectangular  regions  taken 
from  Shriver  et  at.  [2012].  For  each  simulation,  we  compute  the  area-averaged  root-mean-square  amplitude 
(RMSA)  over  the  globe,  in  each  of  the  five  subregions,  and  outside  the  subregions  (see  Table  2).  This  latter 
calculation  provides  a  crude  indication  of  how  far  the  internal  tides  propagate  away  from  their  source 
regions. 

The  simulations  without  a  wave  drag  (e.g.,  E053,  E054,  and  E055;  Figures  4e  and  4f)  generate  highly  ener¬ 
getic  internal  tides  as  seen  in  the  global  RMSA  values  (Table  2;  column  1);  in  addition,  the  internal  tides 
propagate  far  from  their  sources  (Table  2;  last  column).  This  is  also  true  of  each  of  the  subregions  (Table  2) 
and  can  clearly  be  noticed  in  Figure  4.  When  the  wave  drag  is  applied  to  only  the  barotropic  component  of 
the  flow  (e.g.,  Figures  4b  and  4d),  the  resulting  internal  tides  also  have  large  average  amplitudes  (and  travel 
greater  distances)  compared  to  the  simulations  with  the  bottom  wave  drag  (e.g.,  Figures  4a  and  4c) 
although  their  amplitudes  are  not  as  large  as  those  in  simulations  without  a  wave  drag.  By  comparing  the 
twin  simulations  of  bottom  (E072,  E051,  E056,  and  E058)  and  barotropic  (E073,  E052,  E057,  and  E059)  wave 
drag  scenarios  in  Table  2,  we  find  that  the  average  amplitudes  are  consistently  higher  in  the  barotropic 
wave  drag  cases  over  the  globe  and  for  each  subregion.  In  addition,  the  waves  travel  farther  in  the  baro¬ 
tropic  wave  drag  cases  relative  to  the  bottom  drag  cases.  We  note  that  simulation  E073,  which  has  the 
smallest  scale  factor  (0.2)  of  the  barotropic  wave  drag  cases,  gives  values  comparable  to  the  simulations 
without  a  wave  drag,  suggesting  that  the  scale  factor  0.2  is  too  small  to  differentiate  it  from  the  no  drag 
scenarios. 
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Figure  4.  Amplitude  of  the  M2  surface  elevations  from  the  steric  component,  dominated  by  the  internal  tides,  of  the  sea  surface  height  in  HYCOM  simulations  (a)  E051;  with  wave  drag 
(scale  factor  =  0.5)  applied  to  the  bottom  flow,  (b)  E052;  with  wave  drag  (scale  factor  =  0.5)  applied  to  only  the  barotropic  flow,  (c)  E058;  with  wave  drag  (scale  factor  =  1 .0)  applied  to 
the  bottom  flow,  (d)  E059;  with  wave  drag  (scale  factor  =  1 .0)  applied  to  only  the  barotropic  flow,  (e)  E053;  without  wave  drag,  (f)  E055;  without  wave  drag  but  with  quadratic  bottom  fric¬ 
tion  increased  by  about  100  times  along  the  continental  shelves.  The  amplitudes  are  computed  from  3  months  of  SSH  output. 


3.2.2.  Comparison  to  Along-Track  Altimeter  Data 

Figure  5  displays  maps  of  the  SSH  manifestation  of  the  internal  tides  in  HYCOM  and  in  along-track  altimeter 
data.  The  HYCOM  results  shown  in  Figure  5  all  utilize  3  months  of  model  output.  Qualitatively,  we  find  that 
the  default  configuration  with  a  scale  factor  of  1.0  (Figure  5c;  simulation  E058)  exhibits  greater  resemblance 
to  the  altimetric  tides  than  the  other  scenarios  with  respect  to  both  the  scale  of  the  tidal  amplitudes  and 
their  spatial  distribution.  However,  this  should  not  be  taken  as  indicating  that  simulation  E058  yields  the 
most  accurate  baroclinic  tides  because  the  results  depend  upon  the  length  of  time  series  used  as  we  discuss 
shortly.  A  notable  difference  is  that  the  internal  tides  in  the  Atlantic  Ocean  are  somewhat  weaker  in  the 
default  HYCOM  configuration  (Figure  5c)  than  in  the  altimeter  dataset.  Reasons  for  such  discrepancies  have 
been  expounded  upon  in  Shriver  et  al.  [2012]  and  range  from  errors  in  the  barotropic  tidal  flows  in  HYCOM, 


Table  2.  Area-Averaged  Root-Mean-Square  Amplitudes  (in  cm)  of  the  Steric  Sea  Surface  Heights  Over  the  Globe  (Global)  and  in  the  Five 
Subregions  Shown  in  Figure  4  Using  3  Month-Long  SSH  Output 


Simulation 

Global 

Hawaii 

East  of 
Philippines 

Tropical 

South  Pacific 

Tropical 

SW  Pacific 

Madagascar 

Rest  of 
World  Ocean; 

E053 

1.153 

1.967 

2.751 

2.130 

2.707 

2.215 

0.847 

E054 

1.144 

1.963 

2.687 

2.075 

2.677 

2.200 

0.848 

E055 

1.142 

1.977 

2.731 

2.058 

2.660 

2.181 

0.843 

E072 

1.021 

1.819 

2.483 

1.894 

2.388 

2.007 

0.736 

E073 

1.154 

2.024 

2.710 

2.091 

2.738 

2.224 

0.851 

E051 

0.891 

1.658 

2.234 

1.634 

2.016 

1.802 

0.632 

E052 

1.067 

1.923 

2.533 

1.942 

2.524 

2.122 

0.775 

E074 

1.090 

1.975 

2.572 

1.954 

2.558 

2.114 

0.800 

E056 

0.828 

1.548 

2.055 

1.508 

1.869 

1.719 

0.588 

E057 

1.060 

1.950 

2.511 

1.903 

2.453 

2.122 

0.771 

E058 

0.760 

1.465 

1.893 

1.374 

1.649 

1.614 

0.536 

E059 

1.011 

1.902 

2.410 

1.804 

2.343 

2.059 

0.729 

aThe  area-averaged  amplitude  for  the  rest  of  the  world  ocean  outside  the  five  hot  spot  regions. 
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Figure  5.  Amplitude  (cm)  of  M2  internal  tide  in  (a)  along-track  altimeter-based  analyses,  and  in  HYCOM  simulations  (b)  E051;  with  wave  drag  (scale  factor  =  0.5)  applied  to  the  bottom 
flow,  (c)  E058;  with  wave  drag  (scale  factor  =  1 .0)  applied  to  the  bottom  flow,  (d)  E059;  with  wave  drag  (scale  factor  =  1 .0)  applied  to  only  the  barotropic  flow,  (e)  E053;  without  wave 
drag,  (f)  E055;  without  wave  drag  but  with  quadratic  bottom  friction  increased  by  about  1 00  times  along  the  continental  shelves.  The  amplitudes  of  the  HYCOM  simulations  are  com¬ 
puted  from  3  months  of  SSH  output. 


which  generate  internal  tides,  to  mesoscale  contamination  in  the  altimeter  data.  In  the  absence  of  a  wave 
drag  on  the  bottom  flow,  we  find  the  internal  tides  to  be  too  energetic  and  to  travel  too  far  from  their 
source  regions  (compare  Figures  5e  and  5f  to  Figure  5a). 


Wave  Drag  Scale  Factor 


January 

February 

March 

All  3  months 
Altimeter 
Bottom  wave  drag 
Barotropic  wave  drag 


Figure  6  affirms  that  the  scenarios  with  wave  drag  acting  on  the  bottom  flow  yield  smaller  discrepancies 
between  FIYCOM  and  the  along-track  altimeter  data  than  do  their  counterparts  with  drag  acting  on  only 

the  barotropic  flow.  The  RMSA  in  the 
barotropic  wave  drag  simulations  is 
reduced  by  0.077  cm  between  no  drag 
and  drag  scale  factor  1.0,  with  an 
approximately  linear  dependence  on 
the  drag  scale  factor  (Figure  7;  blue 
curves),  presumably  due  to  the  corre¬ 
sponding  decrease  in  barotropic  tide 
signal  strength  (Figure  3)  leading  to 
the  generation  of  weaker  internal  tides. 
Flowever,  RMSA  in  the  bottom  wave 
drag  simulations  is  reduced  by  the 
larger  amount  of  0.24  cm  (Figure  7;  red 
curves)  between  no  drag  and  drag 
scale  factor  1.0,  indicating  that  the  bot¬ 
tom  wave  drag  effect  on  internal  tides 
is  due  both  to  a  reduction  in  the  baro¬ 
tropic  tides  (which  generate  the  baro- 

Figure  6.  Globally  averaged  root-mean-square  amplitude  (RMSA),  of  HYCOM  bar-  clinic  tides)  and  to  a  direct  damping  of 
oclinic  tidal  elevations  from  all  simulations  versus  the  along-track  altimeter  value.  baroclinic  tides  after  they  are  gen- 

The  scale  factor  of  the  wave  drag  for  each  simulation  is  shown  on  the  top  axis  ,  ,  .  . 

and  the  bottom  and  barotropic  drag  simulations  are  denoted  by  red  and  blue  erated.  The  green  Curve  in  Figure  7b 

vertical  dotted  lines,  respectively.  shows  the  difference  in  RMSA  between 
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Wave  Drag  Scale  Factor  Wave  Drag  Scale  Factor 

Figure  7.  (a)  Globally  averaged  root-mean-square  amplitude  (RMSA),  of  HYCOM  baroclinic  tidal  elevations  for  the  barotropic  and  bottom 
drag  simulations  shown  in  Figure  6  (from  3  months  of  analysis),  (b)  The  difference  in  RMSA  between  the  barotropic  and  bottom  drag  simu¬ 
lations  in  Figure  7a  and  their  deviations  from  the  no  drag  simulation. 


the  barotropic  and  bottom  wave  drag  scenarios,  increasing  by  0.17  cm  between  no  drag  and  drag  scale  fac¬ 
tor  1.0.  In  scenarios  without  a  wave  drag  (E053),  even  those  with  enhanced  quadratic  friction  along  the 
shelves  (E054  and  E055),  the  internal  tides  are  too  energetic  and  travel  too  far  from  the  hot  spot  regions 
(see  Figure  6  and  Table  3).  As  can  be  seen  in  Table  3,  the  conclusions  drawn  above  from  the  global  analysis 
also  apply  to  the  five  subregions.  We  analyzed  a  few  1/25°  FIYCOM  simulations  using  1  month  of  SSFH  out¬ 
put,  and  found  that  the  internal  tide  results  are  consistent  with  results  from  our  1/12.5°  simulations.  More 
detailed  results  from  the  1/25°  simulations  will  be  presented  in  the  future;  we  do  not  further  discuss  them 
here  for  the  sake  of  conciseness,  simplicity,  and  computational  expediency. 


Table  3.  Area-Averaged  Amplitudes  (in  cm)  of  M2  Internal  Tides  From  HYCOM  and  Along-Track  Altimetry  (Altimeter)  Over  the  Five 
Subregions  Shown  in  Figure  4a 


Simulation 

Hawaii 

East  of 
Philippines 

Tropical 

South  Pacific 

Tropical 

SW  Pacific 

Madagascar 

Rest  of 
World  Ocean1 

Altimeter 

0.805 

0.837 

0.843 

0.758 

0.715 

0.371 

E053 

1.294 

1.530 

1.270 

1.234 

1.223 

0.516 

E054 

1.291 

1.503 

1.239 

1.205 

1.218 

0.520 

E055 

1.304 

1.547 

1.225 

1.203 

1.210 

0.519 

E072 

1.203 

1.386 

1.144 

1.086 

1.078 

0.443 

E073 

1.329 

1.507 

1.248 

1.252 

1.230 

0.520 

E051 

1.094 

1.227 

0.991 

0.912 

0.952 

0.377 

E052 

1.304 

1.470 

1.179 

1.188 

1.187 

0.493 

E074 

1.302 

1.450 

1.165 

1.186 

1.166 

0.497 

E056 

1.029 

1.104 

0.924 

0.847 

0.904 

0.349 

E057 

1.285 

1.406 

1.136 

1.124 

1.171 

0.477 

E058 

0.971 

1.028 

0.847 

0.748 

0.830 

0.319 

E059 

1.252 

1.349 

1.076 

1.095 

1.135 

0.454 

aThe  HYCOM  results  employ  3  month-long  SSH  output.  All  computations  are  for  depths  deeper  than  1 500  m. 
h'he  area-averaged  amplitude  for  the  rest  of  the  world  ocean  outside  the  five  hot  spot  regions. 


ANSONG  ETAL. 


OPEN  OCEAN  INTERNAL  TIDE  DAMPING 


6067 


10.1 002/201 5JC01 0998 


©AGU  Journal  of  Geophysical  Research:  Oceans 


Figure  6  also  shows  that  though  there 
is  little  monthly  variance  in  the  globally 
averaged  amplitudes,  using  3  month¬ 
long  time  series  for  the  analysis  yields 
values  that  are  reduced  by  about 
9-15%  from  results  computed  from 
1  month-long  SSH  output.  This  sug¬ 
gests  that  some  of  the  discrepancies 
between  the  model  and  along-track 
altimeter  data  may  be  attributed  to 
the  length  of  model  output  used 
in  the  analysis.  To  further  investigate 
the  effect  of  the  length  of  time  series 
on  the  results,  and  to  determine 
which  simulation  yields  the  most 
accurate  internal  tides,  we  run  four 
simulations  out  longer,  to  obtain  a 
total  of  10  months  of  SSH  hourly 
output.  We  run  the  simulations  for 
only  10  months  because  we  only 
have  hourly  NAVGEM  forcing  over 
this  period.  In  this  exercise,  we  con¬ 
sidered  the  no  drag  simulation,  E053, 
and  three  simulations  with  bottom 
wave  drag:  E051  (drag  scale  0.5),  E056  (drag  scale  0.7),  and  E058  (drag  scale  1.0).  Because  the  M2 
internal  tide  amplitudes  in  the  four  simulations  are  still  decreasing  after  10  months  of  analysis  (Figure 
8),  we  also  analyzed  an  older  FIYCOM  simulation  (Expt-18.5;  Figure  8)  from  which  5  year-long  SSH  out¬ 
put  was  saved.  We  find  that  the  internal  tide  amplitudes  in  Expt-18.5  converge  to  an  asymptotic  value 
after  about  4  years,  reducing  from  RMSA  value  of  0.43  cm  at  10  months  to  an  average  value  of  0.41 
cm  after  4  years,  a  reduction  in  amplitude  of  about  5%.  To  obtain  rough  estimates  of  RMSA  from  the 
four  other  simulations  considered  in  Figure  8,  we  use  the  Expt-18.5  results  as  a  correction  to  them, 
suggesting  that  the  four  simulations  E053,  E051,  E056,  and  E058  will  likely  asymptote  to 
0.62,0.48,0.45,  and  0.42  cm,  respectively,  after  about  4  years.  The  estimated  asymptotic  values  for 
these  four  simulations  are  shown  at  month  60  in  Figure  8.  This  exercise  suggests  that  the  along-track 
altimeter  value  of  0.46  cm  will  lie  between  that  of  simulation  E056  (bottom  wave  drag;  scale  factor 
0.7)  and  E051  (bottom  wave  drag;  scale  factor  0.5)  with  E056  yielding  a  value  slightly  closer  to  the 
altimeter  value  than  E051.  Because  the  RMSA  values  of  E056  and  E051  are  both  very  close  to  the 
altimeter  value  and  because  the  corrections  discussed  above  are  based  on  results  from  an  older  simu¬ 
lation,  it  is  still  not  clear  which  of  the  two  simulations  is  superior  to  the  other.  Flowever,  both  simula¬ 
tions  are  more  accurate  than  the  older  simulation,  Expt-18.5,  and  importantly  for  this  paper,  both  are 
more  accurate  than  the  no-wave-drag  simulation  (E053).  Figure  8  shows  that  both  the  record  duration 
and  the  drag  impact  the  coherent  internal  tide  amplitudes.  It  is  not  clear  that  the  impact  of  record 
duration  is  due  to  an  estimation  problem  in  the  presence  of  background  noise  or  to  nonstationarity. 
The  results  on  the  impact  of  record  duration  shown  in  Figure  8  are  consistent  with  those  of  Zhou 
et  al.  [2015],  which  indicate  that  the  internal  tide  variance  from  an  hourly  sampled  mooring  record 
flattens  out  in  about  3  years  (their  Figure  4;  open  circles),  as  well  as  with  related  results  in  Nash  et  al. 
[2012], 

4.  Summary  and  Conclusions 

In  order  to  roughly  constrain  the  damping  of  low-mode  internal  tides,  we  have  conducted  numerical 
simulations  with  a  high-resolution  eddy-resolving  global  ocean  circulation  model  (FIYCOM)  forced  by 
the  M2  tidal  constituent  and  atmospheric  fields.  Various  damping  scenarios  with  different  strengths  of 
parameterized  internal  wave  drag  or  quadratic  bottom  friction  are  performed.  We  ran  simulations  in 


Figure  8.  Globally  averaged  root-mean-square  amplitude  (RMSA),  of  HYCOM  bar- 
oclinic  tidal  elevations  and  along-track  altimeter  value.  The  HYCOM  results  are 
dependent  on  drag  strength  (see  legend)  as  well  as  on  the  length  of  the  time 
series  used.  The  x  axis  is  a  base  2  log  scale.  Expt-18.5  is  an  older  HYCOM  experi¬ 
ment  from  which  several  papers  have  been  written  [e.g.,  Shriver  et  al.,  2012]  and 
from  which  5  year-long  SSH  output  was  stored.  The  amplitudes  shown  at  month 
60  in  our  current  simulations  are  rough  estimates  based  on  the  change  in  Expt- 
18.5. 
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which  a  parameterized  wave  drag  scheme  was  either  applied  to  the  bottom  flow  (referred  to  as  "bot¬ 
tom  wave  drag"  or  the  "default  configuration"  here)  or  to  the  barotropic  component  of  the  flow 
(referred  to  as  "barotropic  wave  drag").  We  argue  that  parameterized  wave  drag  represents  the  break¬ 
ing  of  unresolved  high  modes  by  barotropic  and  low-mode  baroclinic  flows  over  open  ocean  rough 
topography.  In  another  simulation  (designated  here  as  E053),  the  wave  drag  is  completely  turned  off 
and  the  only  drag  on  the  flow  is  quadratic  bottom  friction.  Finally,  to  test  whether  the  propagating 
low-mode  internal  tides  dissipate  mostly  on  the  continental  shelves,  we  turn  off  the  wave  drag  (or 
apply  it  only  to  the  barotropic  flow)  and  at  the  same  time  substantially  increase  the  quadratic  bot¬ 
tom  friction  along  the  shelves.  We  consider  two  cases  (simulation  E054  and  E055)  in  which  the  bot¬ 
tom  friction  in  shelf  regions  is  increased  by  about  25  and  100  times,  respectively,  above  the 
canonical  deep  ocean  value  of  0.0025.  To  test  the  plausibility  of  the  various  damping  scenarios,  the 
resulting  barotropic  tides  in  each  simulation  are  compared  to  the  data-assimilative  TPX08-atlas  solu¬ 
tion.  The  baroclinic  tides  are  compared  to  about  17  years  of  along-track  altimeter  data. 

The  barotropic  tides  in  our  simulations  with  bottom  wave  drag  (especially  simulation  E051)  compare 
reasonably  well  with  TPX08.  Scenarios  without  a  wave  drag  in  the  open  ocean  result  in  barotropic 
tides  that  are  more  energetic  in  HYCOM  than  in  TPX08.  For  the  barotropic  tides,  we  find  that  the  dif¬ 
ferences  in  the  root-mean-square  error  (RMSE)  between  the  simulations  are  relatively  small  (about 
2  cm). 

The  discrepancies  between  the  baroclinic  tides  and  the  altimeter  data  vary  more  significantly  in  the  tested 
scenarios.  In  the  runs  analyzed  for  this  paper,  the  simulations  with  wave  drag  acting  on  the  bottom  flow 
result  in  the  most  accurate  baroclinic  tides.  In  simulations  in  which  a  wave  drag  is  turned  off  entirely,  or  in 
which  it  acts  on  only  the  barotropic  flow,  the  baroclinic  tides  are  overly  energetic  and  propagate  too  far 
from  the  hot  spot  regions,  even  in  the  simulations  (E054  and  E055)  in  which  quadratic  bottom  friction  is 
greatly  increased  in  shelf  regions.  All  of  this  suggests  that,  although  some  internal  tidal  damping  takes  place 
along  the  shelves  [Nash  et  at.,  2004;  Martini  et  at.,  2011;  Kelly  et  al.,  201 3;  Waterhouse  et  at.,  201 4],  substantial 
low-mode  damping  likely  occurs  in  the  deep  ocean.  Waterhouse  et  al.  [2014]  estimate  that  about  30%  of  the 
internal  tide  energy  reaches  into  the  shelves  from  the  open  ocean.  Their  30%  result  is  not  inconsistent  with 
our  conclusion  that  a  substantial  damping  also  takes  place  in  the  open  ocean.  Analysis  of  an  older  FIYCOM 
simulation  (M.  C.  Buijsman  et  al.,  submitted  manuscript,  2015)  indicates  that  bottom  wave  drag  damps 
about  40%  of  internal  tide  energy,  and  that  inferred  residual  numerical  damping  may  dissipate  about  30- 
40%  of  the  internal  tide  energy.  More  analysis,  including  of  the  residual  damping,  needs  to  be  performed 
before  a  final  accounting  can  be  made  of  the  partitioning  of  internal  tide  damping  between  various 
mechanisms. 

We  have  found  that  the  duration  of  model  output,  as  well  as  the  damping  strength,  affects  the  comparison 
of  modeled  internal  tides  to  satellite  altimetry.  Figure  8  demonstrates  that  the  globally  averaged  internal 
tide  amplitudes  in  our  simulations  require  about  3-4  years  to  converge.  The  long  time  required  for  conver¬ 
gence  is  consistent  with  results  in  Nash  et  al.  [2012]  and  Zhou  et  al.  [2015],  who  both  analyzed  tide  records 
at  individual  locations.  It  is  not  clear  whether  this  is  due  to  an  estimation  problem  in  the  presence  of  back¬ 
ground  noise  or  to  nonstationarity,  and  additional  studies  are  required  to  determine  the  underlying  cause. 
We  note  that  an  increasing  number  of  internal  tide  studies  conclude  that  internal  tide  incoherence  is  a  sub¬ 
stantial  effect.  Zaron  [2015],  for  example,  finds  that  about  one-third  of  the  global  semidiurnal  internal  tide 
sea  surface  height  variance  is  incoherent.  Using  5  year-long  SSFH  output  from  an  older  FIYCOM  experiment, 
we  estimate  asymptotic  values  of  the  M2  amplitudes  for  four  of  our  simulations  (E053,  E051,  E056,  and  E058; 
Figure  8),  and  find  that  simulations  E051  (bottom  wave  drag  with  scale  factor  0.5)  and  E056  (bottom  wave 
drag  with  scale  factor  0.7)  yield  the  most  accurate  baroclinic  tides. 

We  reiterate  that  our  simulations  represent  various  limiting  cases  of  the  aforementioned  damping  mecha¬ 
nisms.  Other  plausible  scenarios  for  damping  low-mode  internal  tides,  such  as  upper-ocean  wave-wave 
interactions,  have  not  been  tested  in  this  study  and  will  instead  be  investigated  separately.  It  remains  to  be 
seen  whether  upper-ocean  wave-wave  interactions  represent  an  open  ocean  damping  mechanism  that  is 
equally  plausible  as  bottom  wave  drag.  In  the  actual  ocean,  of  course,  more  than  one  mechanism  is  likely  to 
be  important.  We  leave  it  to  future  work  to  examine  model  runs  having  more  than  one  open  ocean  damp¬ 
ing  mechanism  acting  simultaneously,  an  endeavor  considered  to  be  too  complicated  for  the  present 
paper. 
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In  future  work,  we  will  examine  in  more  detail  the  impacts  of  record  duration  and  model  resolution.  We  will 
also  examine  the  energy  budgets  in  simulations  with  different  drag  strengths. 
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